Quantum channels in random spin chains 
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We study the entanglement between pairs of qubits in a random antiferromagnetic spin- 1/2 chain at zero 
temperature. We show that some very distant pairs of qubits are highly entangled, being almost pure Bell states. 
Furthermore, the probability to obtain such spin pairs is proportional to the chain disorder strength and inversely 
proportional to the square of their separation. 
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I. INTRODUCTION 

The concept of entanglement between two systems dates 
back to the birth of quantum mechanics [0, 0. In a certain 
sense, entanglement can be associated with the counterintu- 
itive nonlocal behavior of two quantum systems: namely, the 
appearance of nonlocal correlations between them even af- 
ter their mutual interaction is switched off. Such a feature 
has puzzled our minds, and explorations of that concept have 
opened new fields of research such as quantum computation 
and quantum information |0] . In particular, many useful quan- 
tum communication protocols such as superdense coding 01, 
teleportation and quantum state sharing 10,0] require en- 
tangled states (quantum channels) as a key ingredient for their 
successful implementation. 

In view of that, the search for physical systems where en- 
tangled states can be found turned to be of utmost importance 
not only for academic reasons, but also for practical imple- 
mentation of the new technology arising from these new fields 
of research. Promising physical systems in which entangle- 
ment may be abundant are quantum spin chains. For a long 
time one has intensively studied such systems mainly because 
of their interesting magnetic properties. Therefore, a great 
amount of knowledge about the physics of such many-body 
systems was obtained, paving the way to an investigation of 
the entanglement properties of spin chains H). A fruitful re- 
sult of this investigation, for example, is the discovery of the 
role of entanglement in quantum phase transitions. At the 
critic alpoint, it is known that the bipartite entanglement en- 
tropy H of the ground state diverges logarithmically with the 
size of the subsystem. In addition to this, it was shown that 
the bipartite entanglement is proportional to a quantity called 
central charge, which dictates the conformal field theory uni- 
versality class of the chain lioll . 

Almost all previous studies dealt with ordered spin chains. 
Recently, nevertheless, it was shown that the bipartite entan- 
glement of a random antiferromagnetic (AF) Heisenberg spin- 
1/2 chain al so g rows logarithmically with the size of the sub- 
system ifTTl Il2ll . A similar result was also derived for a ran- 



dom transverse-field Ising chain at the critical point. In both 
cases the 'effective' central charge is less than the one calcu- 
lated for their clean and ordered counterparts ifTlll . In other 
words, bipartite entanglement entropy is decreased when dis- 
order is added into the system. In view of such a result, one 
may naively think that disorder is an undesired feature in spin 
chains when one is concerned with implementation of any 
quantum information protocol. Here we show, interestingly, 
that disorder provides a ground state with remarkable entan- 
glement properties which can be properly harnessed to create 
almost maximally entangled Bell states, the ultimate quantum 
channel needed in many quantum communication protocols. 

In this work we study the entanglement properties of long- 
distance pairs of qubits in the random AF spin- 1/2 chain. For 
that purpose we determine for some suitable pairs in the chain 
how close they are to a Bell state; i.e., we compute the fi- 
delity of those pairs with a singlet. In addition, we exhaus- 
tively characterize their entanglement properties computing 
their negativity fl3ll . logarithmic negativity lfl3ll . and concur- 
rence (l4j|. Moreover, we relate all these quantities to each 
other. We show that some rare pairs display strikingly high 
fidelities, qualifying them as good quantum channels. Finally, 
we show that the fidelity distribution of such rare pairs does 
not depend on the size of the chain and is proportional only to 
the disorder strength. 

The remainder of this paper is organized as follows: We 
introduce the model, compute the fidelity, and relate it to the 
negativity, logarithmic negativity, and concurrence in Sec. ITU 
In Sec. [HI] we present our numerical data, quantifying the in- 
teresting pairwise entanglement properties of the disordered 
system, and compare it with the ordered and clean ones. Fi- 
nally, we conclude and discuss the relevance of our results in 
SecHV] 



II. THE MODEL AND ITS PAIRWISE ENTANGLEMENT 
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In this section we introduce the model and compute its pair- 
wise entanglement as a function of the spin-spin correlation 
functions. 
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A. The model 

Our starting point is the random XXZ AF spin- 1/2 chain 
whose Hamiltonian reads 



(2.1) 



where 5; is the usual spin- 1/2 operator related to site i and 
/, > and A, are independent random variables. In all the 
following calculations we use A, = 0, V;. However, we first 
argue that our results are valid for all models satisfying the 
condition — 1/2 < A,- < 1. It was shown that weak anisotropy 
(— 1/2 < A,- < 1) renormalizes to zero in the low-energy limit, 
being thus irrelevant Il5lll6ll . Furthermore, the isotropic point 
A, = 1, V/, is governed by the same infinite-randomness fixed 
point (IRFP) that governs the point A, = 0; i.e., those sys- 
tems quantitatively share the same long-length ground state 
properties IT^l . Thus, from now on it is sufficient to con- 
sider only the case where A, = 0. This choice allows us to 
map the Hamiltonian (12. U to a free fermion one through a 
Jordan-Wigner transformation. Then, we can exactly diago- 
nalize large chains using standard subroutines lfl7ll . 

The phase diagram of the clean counterpart of this system 
(/, = / and A, = A) is well known [18]. For -1 < A < 1 the 
system is a critical Tomonaga-Luttinger spin liquid lfl9ll . For 
A > 1, a spin gap opens and the system belongs to the Ising 
antiferromagnet universality class, and for A < — 1 the system 
is in the gapless Ising ferromagnet universality class. In the 
critical regime, the spin-spin correlation function decays as a 
power law with well-known exponents B20I1 and conjectured- 
exact amplitudes ||2lll . In addition, it was shown that the pair- 
wise entanglement between nearest-neighbor spins is maxi- 
mum at the isotropic point I22I l23ll . However, the pairwise 
entanglement between further neighbors has a maximum near 
A = — 1 i23ll . Moreover, the ground-state pairwise entangle- 
ment is not affected by spontaneous symmetry breaking 12311 . 

In contrast to the clean model, the low-energy physical be- 
havior of all the chains parametrized by Hamiltonian d2.lt in 
which — 1/2 < A, < 1 is completely different, no matter how 
weak the disorder strength is. In the renormalization-group 
(RG) sense, it is said that the chains are governed by the so- 
called IRFP |[l6tl . Within this framework, their ground states 
are a collection of spin singlets whose spin pairs can be arbi- 
trarily faraway from each other [see Fig.[TJa)]. 




Figure 1: (Color online) (a) Pictorial view of the random singlet 
ground state. Bonds connecting spins denote singlet pairs, (b) The 
ground state from the point of view of the highlighted spin pair. 

It is also known that the chance of finding a very distant 



spin singlet pair is rare — viz., ~ \i — j\~ . However, these 
rare spin pairs contribute enormously to the average spin- 
spin correlation function, since they develop correlations of 
order unity. Hence, the mean correlation function decays al- 
gebraically as \i — j'p 2 flrjll . In contrast, a typical faraway 
spin pair is not in a singlet state, developing weak correlation 
~ exp{ — — j\}, meaning a typical correlation function de- 
caying as exp{-^|/-;'|} HI]. 

How close those rare faraway singlets are to a perfect sin- 
glet (or Bell) state is the question we address in this work. 
Moreover, we wish to quantify their entanglement. As we 
show in Sec. [Ill] there are rare and extreme faraway spin pairs 
displaying enormous fidelity and entanglement. This is a re- 
markable property absent in the clean and ordered systems: 
Two spins separated by a distance of order of the chain length 
and interacting indirectly through the spins between them are 
strongly entangled and, therefore, weakly entangled with the 
rest of the chain. 

We now explain how those rare spin pairs are located in 
the chain. Given a disordered chain they can be easily found 
using the decimation rules predicted by the strong-disorder 
renormalization group method iflrj l24ll . One has to search 
for the greatest coupling constant in the chain — say, J2 — 
and then decimate spins 5i and S3 by coupling them in a 
singlet state. The neighbors spins Si and S4 now interact 
with each other through a new renormalized coupling con- 
stant /= JiJi/ [(1 + A2)/2]- The local anisotropic parameter 
also renormalizes to A = A1A3 (1 + A2) /2. Iterating this pro- 
cedure one can find the position of all spin pairs in the chain. 
Since all the coupling constants (including the renormalized 
ones) are random independent variables, the singlet pairs are 
distributed in a random fashion, and thus the system is said to 
be in a random singlet phase. 



B. Entanglement properties 

Having localized the spin pairs we now want to study how 
close they are to a pure singlet — i.e., the Bell state I 1 ? - ) = 
(| +— ) — I — h}) /V2, with |+) and |— ) meaning spin up and 
down, respectively. For this purpose we calculate the fidelity 
Fij between the actual state describing the random singlets py 
and I^P - ), which is given as 



F ij = (m-\ 9ij \m-), 



(2.2) 



where py = tryj(p) is the usual reduced density matrix of spins 
Si and Sj (obtained tracing out all spins but 5, and Sj) and 
p = \g) (g\ is the zero temperature density matrix in which 
|(?) is the ground state of the Hamiltonian (12.11 ). Note that 
Ffj = 1 if py is a perfect singlet and that fy = whenever it 
is orthogonal to I 1 ? - ). 

Using symmetry arguments we can show that py = 

4£ a CffSfSJ + 1 /4, with a = x, y, z and Cff = {§ SfS a j Q ^ 
being the ground state spin-spin correlation. In addition, 
Cfj = Cy for the XXZ AF spin- 1/2 chain. Expanding py 

in the Bell basis, |4>±) = (| + +) ± | - -))/\/2 and pP*) = 



3 



(l + -)±|- 

Py 

where 



-f))/V5, we get 

i+C^(|4> + )(4> + | + |4>-)<<I>-|), (2.3) 



■ qj. 



(2.4) 



However, we are also interested in a complete quantifica- 
tion of the entanglement between spins 5, and Sj. There- 
fore, we compute its negativity lfl3ll which can be defined as 
9&j = I IP/jl I — 1> where | |A|| = ti^VA'A) is the trace norm of 
operator A and 7] means partial transposition, in the sense that 
the transposition operation is performed only in spin Sj l25ll . 
The normalization for 5Y} ; - was chosen such that it is equal to 
1 for a maximally entangled state (e.g., Bell state). A direct 
calculation using Eq. ( 12.3b gives 



Ha 



0. 



iffy < 1/2, 



2Fij-l, if Ffj > 1/2. 



(2.5) 



We remark that whenever Fy > 1/2 the two spins are entan- 
gled and that the greater the fidelity the greater the entangle- 
ment. The entanglement threshold Fy = 1/2 will be an impor- 
tant quantity in the discussions below. Note that the negativity 
is a function of the fidelity, implying that knowledge of the lat- 
ter allows us to obtain the negativity through Eq. ( 12.5b . 

Moreover, a complete characterization of the entanglement 
properties is desirable. Interestingly, for the models described 
by the Hamiltonian d2.lt . we find that the logarithmic nega- 
tivity £y lfl3ll and the concurrence Cy ifTill can be expressed 
in terms of the fidelity: 



£y = log 2 



0, 

log 2 (2fy) 



if F u < 1/2, 
iffy > 1/2, 



and 



Qj = max{0, A-i - 14- A, 3 - A4} = !SQj, 



(2.6) 



(2.7) 



where the X's are the square roots of the eigenvalues of 
the matrix py 



&-®Gjj p*j {Oj ® OjJ arranged in decreas- 
ing order (26J, o y is the usual Pauli matrix, and p* is 
the complex conjugate of p. The logarithmic negativity is 
an additive entanglement measure for mixed states not re- 
covering the von Neumann entropy for pure states in gen- 
eral lll3ll However, the entanglement of formation iflill . £/ = 

-x\og 2 x- (1 -jt)log 2 (1 -x), withx = 1/2 + Vl - C 2 /2, is 
a monotonous function of the concurrence and does reduce to 
the von Neumann entropy for pure states. 

We thus conclude that, due to the symmetries of the Hamil- 
tonian d2.lt . knowledge of the fidelity is sufficient to com- 
pletely quantify the entanglement between any pair of qubits 
in the chain. 



in. NUMERICAL RESULTS 

We now report our numerical study on the fidelity between 
spin pairs. Here we employed chains with coupling constants 
distributed according to 



P(J) 



(1 



a 



ft, 



l-a 



■r 



(3.1) 



where < / < £2o and the disorder parameter a < 1 . Hence, 
greater a means greater disorder and Ho is the energy 
scale l27tl . We have considered three different disorder pa- 
rameters a = 0, a = 0.3, and a = 0.6 and chains of length 
Lo = 100, 200, 400, and 800 with periodic boundary condi- 
tions. 

In order to emphasize the difference between the clean and 
disordered systems, we plot in Fig. [2] the histogram of fideli- 
ties for all spin pairs separated by odd lengths greater than a 
cutoff l c = Lq/6 1 28]. As shown in panel (a), the fidelities of 
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Figure 2: (Color online) Normalized histogram of fidelities for those 
spin pairs separated by distances greater than a cutoff l c = Lo/6 for 
the disordered [a = 0.3, leftmost solid (black) and dashed (green) 
curves] and the clean and ordered chains [rightmost solid (orange) 
and dashed (purple) curves], and for length chains Lo = 100 (solid 
lines) and Lo = 200 (dashed lines), (a) The typical value of F is 
clearly greater in the clean chain than in the disordered one. (b) 
However, as revealed by the logarithmic scale, there are some rare 
spin pairs (rare events) displaying high fidelities in the disordered 
chain, while in the clean one the fidelities are sharply constrained to 
values less than 1/2. A total of 6000 different disordered chains were 
used to build each one of these histograms. 

those spins in the disordered chain are typically less than those 
in the ordered one. This is because in the disordered system 
the typical correlation function decays ~ exp{— j\}; 
hence, the typical value of Fy approaches 1 /4 rapidly when Lo 
is increased, while in the clean system the transversal and lon- 
gitudinal correlation function decays as |z — with r| = 1/2 
and 2, respectively Ir29ll . Moreover, as Cf" depends only on 
— in the ordered system, the fidelities are thus sharply 
bounded as dictated by Eq. d2.4t and clearly shown in panel 
(b). As also shown in panel (b), on the other hand, the fi- 
delities in the disordered chain are not sharply bounded from 
above. It is due to the fact that C"" is distributed for a fixed 
\i — j\. It is worth mentioning that there exist some rare spin 
pairs displaying strikingly high fidelities. Although they are 
very rare, their fidelities are above the threshold 1/2 below 
which there is no entanglement [see Eqs. d2.5b -( r2~7| i1. In con- 
trast, there is no distant spin pair displaying such high fidelity 
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in the clean chain. 

We now focus our attention on those rare and distant pairs. 
In order to characterize their entanglement, we choose one 
spin pair per chain to represent them. There are two natural 
choices: (i) the longest and (ii) the last decimated spin pairs. 
The first choice is interesting because we are focusing on pairs 
whose spins are extremely distant. However, we choose the 
second one by the following reasons: (i) it is easier and faster 
to implement numerically, (ii) there is a great chance the last 
decimated pair be the longest one, and (iii) as we shall see 
below, we focus our attention only on pairs separated by dis- 
tances of order of the chain length. We thus proceed as fol- 
lows. 

We entirely decimate a chain following the RG rules l24ll in 
order to find the position of the last spin pair. Then, we com- 
pute their spin-spin correlations exactly [17]. With this infor- 
mation we calculate the fidelity between this last pair and a 
singlet as given in Eq. ( 12.41 ). In Fig. [3] we plot the fidelity dis- 
tribution Q (F) for chains with a = 0, 0.3, and 0.6 [see pan- 
els (a), (b), and (c), respectively] for different chain lengths. 
For each length and disorder parameter, we computed a total 
of 10000 fidelities in order to build the distributions Q(F). 
In panels (a)-(c) we retained only those pairs separated by 
a distance between Lq/6 and Lo/2 (around 70% of the last 
spin pairs follow this requirement, independent of disorder 
parameter and chain length). In all the cases studied, there 
is no difference between the histograms for different chain 
lengths. We thus conclude that the probability of finding/ar- 
away and highly entangled spin pairs does not depend on the 
chain length and increases with the disorder strength. 
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Figure 3: (Color online) Fidelity distribution Q(F) for the chains 
with disorder parameter (a) a = 0, (b) a = 0.3, and (c) a = 0.6, con- 
cerning the last decimated spin pairs separated by distances greater 
than l c = Lo /6. The same in panel (d) for a fixed length Lq = 200 
and different disorder parameters. Open (solid) symbols refer to spin 
pairs separated by distances greater (shorter) than l c . The vertical 
dashed line highlights the threshold F = 1/2 below which there is no 
entanglement. 

Moreover, looking at panels (a)-(c) and especially panel (d) 



of Fig. [3] we clearly see that an increasing disorder increases 
the probability of obtaining a spin pair with greater fidelity 
and greater entanglement. We also highlight that for greater 
disorder [panel (c)] we have a reasonable probability of find- 
ing pairs of spins with fy « 1. Those pairs can be considered, 
for all practical purposes, perfect singlets. This is in contrast 
with ordered and clean chains in which pairwise entanglement 
is found only between nearest neighbors II 8l l23l 13011 . 

We want now to draw some attention to the following 
caveat. The mean value of the fidelities concerning all spin 
pairs separated by distances of order of the chain length ap- 
proaches 1/4 as the chain length increases. However, there 
are always some rare spin pairs displaying enormous entan- 
glement [as shown in Fig. |2fb)]. Moreover, the probability 
of finding them is ~ \i — j\~ 2 ; i.e., it depends only on their 
separation lfl6ll . Precisely, for \i — j\ 3> 1, the probabity of 
finding them equals 2/(3 \i — j\ 2 ) l3~lll . Furthermore, the fi- 
delity distribution concerning only those rare spin pairs Q (F) 
is completely different. It does not depend on the separation 
between them, and its mean value increases with the disorder 
strength [as shown in Figs.[3ja)-[2c)]. 

Our results are in agreement with the localized nature of the 
ground state. Spin pairs separated by distances greater than 
the corresponding order-disorder crossover length shall share 
correlations mediated by virtual polarizations of the spin sin- 
glets between them 13211 . Note that the fidelity distributions 
of those last pairs separated by distances shorter than Lq/6 is 
slightly shifted to the right [see Fig.[3ld)]; i.e., the mean value 
of the fidelity is greater. This result does not mean there is a 
length dependence on Q (F) since there are many other short- 
length spin pairs earlier decimated without such high fidelity. 
Moreover, there are pairs separated by distances less than the 
ordered-disordered crossover length, meaning we cannot sim- 
ply compare such distributions. 

Besides the existence of a reasonable number of almost 
pure singlets for strong disorder, we also have a remarkable 
ensemble of states suited for the implementation of distilla- 
tion protocols lO . A distillation protocol can be seen as a 
local procedure by which we can obtain M sin glet s from an 
ensemble ofN>M mixed states. One can show IB3I1 that only 
mixed states py with a fidelity, as given by Eq. d2.2l ). greater 
than 1/2 can be employed for distillation purposes. And this 
is exactly the scenario we have in random AF spin- 1/2 chains. 
For example, for a disorder parameter a = 0.3 we have almost 
90% of the last decimated singlets with fy > 1/2 and we can 
find many pairs of spins more than 100 sites apart with fideli- 
ties above this threshold. For a = 0.6 we achieve rates greater 
than 97%. Remarkably, those spins are extremely distant from 
each other — namely, fractions of the chain size. 

Typically, the effective coupling J e ff between distant spins 

is exponentially weak — viz., £2o ex p{~Y\/l' — 71} > where y 
is a constant proportional to the disorder strength but of or- 
der of unit [16]. Analyzing the distribution of J e ff we find 
y« 1.2, 2.0, and 3.5 for a = 0, 0.3, and 0.6, respectively. 
Thus, the typical value of J e /f with respect to a chain with 
disorder parameter a = 0.6 for a pair separated by a distance, 
say, =40 is of order 10~ 10 i2o. With this result, one 

may be concerned about effects of local perturbations on the 
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singlets and also effects of finite temperature. How ever, a l- 
though the typical value of J e tt decreases ~ exp{— y[?--7|}, 
the width of the distribution of In (/<,//) increases ~ — j\, 
meaning there is a finite probability to find a high-fidelity spin 
singlet with relatively strong coupling. Among our data re- 
garding chains with lengths up to 100 sites and a = 0.6, we 
found 5 chains where the distances between the spins were 
greater than 40 sites, the fidelities were greater than 0.75, and 
their effective couplings were greater than 10~ 6 £2o. From the 
energetic point of view, less disordered chains are favorable. 
For the chain with a = and same length, we found 1 1 of 
those spin pairs but their effective coupling were greater than 
10- 4 i2 . 

IV. FURTHER DISCUSSIONS AND CONCLUSION 

We now put our results in a more general context. Consider 
first the random AF dimerized spin chains whose Hamiltonian 
is given by Eq. ( 12. U with different coupling constant distri- 
butions between odd and even bonds. In the Griffiths phase 
near criticality (i.e., the even and odd coupling constant dis- 
tributions are not identical but have great overlap), the strong- 
disorder RG scheme also applies and spin pairs couple in sin- 
glet states l34ll . Thus, random spin pairs resembling perfect 
singlets are also expected. However, the probability of finding 
distant spins in a singlet state is lower when compared to the 
model studied here. Also, it is known that the random AF two- 
leg and zigzag spin- 1/2 ladders renormalize to random critical 
or dimerized spin chains in certain conditions j35ll . This sug- 
gest that the results here presented also apply to those systems. 

Let us now raise an interesting question: if the singletlike 
spin pairs are weakly perturbed, how do they behave after a 
projective measurement? It is reasonable to expect they will 
oscillate coherently with a period inversely proportional to 
their effective coupling constant. This means that spins \i — j\ 
sites away will slowly oscillate with period ~ exp {y]7— /|) 
while the rest of the chain stays unperturbed IBol . Note that 
in this strong-disorder scenario there is no room for the char- 
acteristic "decoherence wave" of the clean system l37ll . The 
mechanism behind such "protection" against decoherence is 
certainly related with the localized aspect of the system, where 
the effective coupling between the spin pair and the rest of the 
chain is weak [see Fig. [Tib)]. 

In conclusion, we have shown that the ground state of ran- 
dom AF spin- 1/2 chains displays remarkable entanglement 
properties. First, it is possible to find long-distance spin pairs 
that are almost perfect singlets. Actually, the more disorder 



the system possesses, the greater is the probability of find- 
ing such pairs, meaning that they become weakly perturbed 
singletlike states. These results are in agreement with what 
is expected from the localized nature of the ground state as 
given by the strong-disorder RG scenario of the random sin- 
glet phase. However, the fidelities and negativities here cal- 
culated provided the first suitable and exact quantification of 
how much such pairs resembles a singlet state. Moreover, 
we have shown that the distribution of the fidelity and en- 
tanglement between distant singlet-like pairs does not depend 
on the chain length, depending only on the disorder. Conse- 
quently, the probability to find a high-fidelity spin pair is equal 
to 2Q(F)/(3 \i — j\ 2 ), provided the distance \i — j\ between 
them be much longer than the corresponding order-disorder 
crossover length. Notice the only dependence on disorder is 
given in Q (F) as shown in Figs. Sa)-[2c). In addition to 
showing that some spin pairs are very good quantum channels 
(fidelity with a singlet wl)we have shown that all the weakly 
perturbed spin pairs can be used to create perfect singlets. For 
strong disorder we have shown that more than 90% of these 
pairs possess fidelities greater than the lower bound 1/2 above 
which quantum distillation protocols, procedures in which a 
few perfect singlets can be obtained from a large ensemble of 
mixed states, become feasible. 

We also point out that the use of local defects in related 
models in order to achieve strongly entangled states was re- 
ported in Ref. ll38ll . However, the entanglement there were 
obtained between near spins (qubits) in contrast to those here 
reported. 

Finally, it has been recently achieved an experimental real- 
ization of the AF spin-5/2 Heisenberg model [39]. The aston- 
ishing experimental precision achieved by the scanning tun- 
neling microscope leads to an unprecedent experimental pos- 
sibility of confronting the theoretical predictions of the effects 
of disorder on low-dimensional magnets. In particular, this 
experiment suggests that in the near future the results here 
presented can be tested and applied to the creation of highly 
entangled quantum channels. 
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